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1.0  EXECUTIVE  SUMMARY 


Under  suitable  eonditions,  monitored  natural  attenuation  (MNA)  ean  be  a  cost-effeetive  strategy 
for  restoring  eontaminated  aquifer  systems  either  as  a  stand-alone  teehnology  or  in  eombination 
with  other  engineered  remedial  aetions.  However,  the  Environmental  Proteetion  Ageney  (EPA) 
guidance  specifically  requires  MNA  to  achieve  site-specific  cleanup  objectives  within  a 
reasonable  time  frame.  Thus,  it  is  necessary  to  provide  estimates  of  cleanup  times  whenever 
MNA  is  proposed  as  part  of  a  cleanup  strategy.  In  response,  a  screening  tool  NAS  was  developed 
for  estimating  time  of  remediation  (TOR)  for  MNA  with  varying  degrees  of  source  area 
remediation.  NAS  is  designed  to  make  complex  analytical  and  numerical  solutions  of  the  TOR 
problem  accessible  to  remedial  project  managers  (RPM)  and  their  contractors  using  site-specific 
remediation  objectives.  Conventional  screening  tools  for  MNA  are  not  designed  to  address 
source  zone  remediation  options  or  simulation  plume  reduction. 

The  Natural  Attenuation  Software  tool  was  co-developed  by  the  U.S.  Navy,  USGS,  and  Virginia 
Tech.  NAS  consists  of  a  combination  of  computational  tools  implemented  in  three  main 
interactive  modules  to  provide  estimates  for:  (1)  target  source  concentration  required  for  a 
plume  extent  to  contract  to  regulatory  limits,  (2)  time  required  for  contaminants  in  the  source 
area  to  attenuate  to  a  predetermined  target  source  concentration,  and  (3)  time  required  for  a 
plume  extent  to  contract  to  regulatory  limits  after  source  reduction.  Natural  attenuation  processes 
that  NAS  models  include  are  advection,  dispersion,  sorption,  non-aqueous  phase  liquid  (NAPE) 
dissolution,  and  biodegradation  of  petroleum  hydrocarbons,  chlorinated  solvents,  or  any  user- 
specified  contaminants  or  mixtures. 

The  objective  of  this  demonstration  is  to  evaluate  the  NAS  capability  to  provide  reasonable 
estimates  of  MNA  cleanup  time  frames  in  a  variety  of  environments  and  sites  throughout  the 
United  States.  The  tool  is  evaluated  by  using  data  from  eight  sites  with  long-term  monitoring 
data  that  encompass  diverse  geologic  and  hydrogeochemical  environments  and  different 
remediation  options.  The  eight  demonstration  sites  are  located  at  Seneca  Army  Depot  and  a 
USGS  study  site  in  New  York,  Naval  Air  Engineering  Station  (NAES)  Eakehurst  in  New  Jersey, 
Hill  Air  Eorce  Base  (AFB)  in  Utah,  Naval  Submarine  Base  (NSB)  Kings  Bay  in  Georgia,  Naval 
Air  Station  Cecil  Eield,  Naval  Air  Station  Pensacola  in  Elorida  and  a  USGS  study  site  in  Alaska. 

Remedial  actions  at  Comprehensive  Environmental  Response,  Compensation,  and  Eiability  Act 
(CERCEA)  sites  must  be  protective  of  human  health  and  the  environment  and  comply  with 
applicable  or  relevant  and  appropriate  requirements  (ARAR).  Drinking  water  standards  provide 
relevant  and  appropriate  cleanup  levels  for  ground  waters  that  are  a  current  or  potential  source  of 
drinking  water.  Drinking  water  standards  include  federal  maximum  contaminant  levels  (MCE) 
and/or  non-zero  maximum  contaminant  level  goals  (MCEG)  established  under  the  Safe  Drinking 
Water  Act  (SDWA),  or  more  stringent  state  drinking  water  standards.  Selection  of  MNA  as  a 
remedy  for  cleanup  of  hazardous  substances  requires  that  estimates  of  time  of  cleanup  be  made 
to  assess  its  feasibility  as  a  remedy. 

NAS  was  effective  in  predicting  the  time  of  stabilization  of  concentrations  at  downgradient 
monitoring  wells  located  relatively  close  to  the  source  (within  100  to  700  ft)  following  source 
remediation  and  a  reduction  in  groundwater  contaminant  concentrations  in  the  source  zone. 
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Accuracy  of  the  solution  relative  to  the  observations  was  impacted  by  dynamic  behavior  of  the 
plume  over  time.  NAS  was  also  effective  in  capturing  concentration  time  trends  of  natural  source 
depletion  of  a  multieomponent  NAPL,  providing  a  prediction  that  was  superior  to  a 
comprehensive  numerical  model  that  was  not  based  on  source  zone  mass  balance. 

One  finding  of  this  demonstration  was  that  NAS  proved  to  be  applieable  to  all  eight  sites, 
independent  of  hydrogeology,  contaminants,  characteristics  of  the  source  zone,  or  engineered 
remedial  action  (ERA).  Therefore,  the  simplifying  assumptions  associated  with  the  analytical 
solutions  and  the  numerieal  source  zone  model  do  not  appear  to  render  NAS  ineffective  but,  in 
fact,  demonstrate  the  applieability  and  utility  of  NAS  to  a  wide  range  of  contaminated  sites.  In 
contrast,  comprehensive  three-dimensional  numerieal  models  constructed  to  simulate  the 
complexities  of  a  groundwater  system  and  features  of  a  plume  subject  to  limited  data  may 
include  unrealistic  boundary  conditions  that  do  not  honor  the  actual  field  conditions. 

An  example  of  an  approach  for  combining  NAS  with  cost  to  completion  analysis  is  presented. 
An  estimate  of  the  life-eycle  costs  of  monitoring  associated  with  MNA  and  natural  souree  zone 
depletion  at  NAS  Ceeil  Field,  Site  3  was  developed  in  conjunction  with  NAS  TOR  estimates. 
Using  the  estimated  TOR  for  triehloroethene  (TCE)  (43  years)  and  naphthalene  (>69  years),  the 
net  present  value  for  monitoring  over  this  period  of  time  is  at  least  $2,333,029.  In  this  case, 
additional  optimization  at  Site  3  was  reeommended  due  to  the  relatively  high  life-cyele  costs  as 
well  as  the  loss  of  efficiency  associated  with  MNA  in  the  source  area  resulting  from  air  sparging. 

An  estimate  of  the  costs  to  implement  NAS  relative  to  a  comprehensive  numerical  solute 
transport  model  was  presented.  The  analysis  demonstrated  that  a  cost  savings  between  $94,650 
and  $126,030  ean  be  realized  by  using  the  NAS  modeling  approaeh.  The  cost  savings  associated 
with  NAS  refieets  a  savings  in  labor  by  a  factor  of  5-6  and  a  reduction  in  data  requirements 
relative  to  the  comprehensive  approach.  If  this  same  cost  savings  was  achieved  at  sites  (204 
total)  on  an  ineremental  basis  over  5  years,  then  the  cumulative  5-year  eost  savings  to  the  DoD 
would  be  $20,900,000.  This  analysis  does  not  factor  in  the  additional  potential  eost  savings  to 
DoD  at  these  sites  by  implementing  an  MNA-based  strategy,  with  or  without  source  zone 
reduction. 

A  methodology  and  tool  for  estimating  the  time  of  remediation  associated  with  MNA  allows 
stakeholders  to  make  informed  decisions  regarding  its  application.  In  addition,  budget 
requirements  for  long-term  monitoring  programs  can  be  forecast  based  on  estimates  of  time 
frames.  This  allows  better  program  planning  to  meet  the  future  needs  of  cleanup  programs,  and 
can  afford  RPM  the  ability  to  conduet  cost  benefit  analyses  when  comparing  source  removal 
with  MNA  options  to  MNA-only  strategies. 
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2.0  TECHNOLOGY  DESCRIPTION 


2,1  TECHNOLOGY  DESCRIPTION  AND  APPLICATION 

MNA  is  the  considered  use  of  naturally  occurring  contaminant  degradation/dispersion/ 
immobilization  processes  to  reach  site-specific  remediation  goals  (EPA,  1996;  1998).  In  current 
engineering  practice,  the  effectiveness  of  MNA  is  evaluated  on  a  site-by-site  basis  by 
considering  three  lines  of  evidence:  (1)  historical  monitoring  data  showing  decreasing 
concentrations  and/or  contaminant  mass  over  time,  (2)  geochemical  data  showing  that  site 
conditions  favor  contaminant  transformation  or  immobilization,  or  (3)  site-specific  laboratory 
studies  documenting  ongoing  biodegradation  processes  (EPA,  1998).  A  variety  of  field  and 
laboratory  methods  for  assessing  these  three  lines  of  evidence  have  been  developed  and  are 
currently  in  use  (Wiedemeier  et  ah,  1999). 

In  concept,  estimating  the  length  of  time  required  for  natural  processes  to  remove  a  particular 
contaminant  from  a  groundwater  system  is  a  mass  balance  problem.  It  requires  estimates  of  the 
rate  of  mass  flux  from  and  depletion  over  time  of  NAPE  or  diffused/sorbed  contaminant  sources. 
It  also  requires  estimates  of  the  rates  of  natural  attenuation  processes  (advection,  hydrodynamic 
dispersion,  biodegradation,  sorption,  plant  uptake,  etc.).  In  practice,  a  deterministic  approach  can 
be  employed  to  quantify  the  mass  balance  of  the  aqueous  phase  plume  and  the  contaminant 
source. 


Widdowson  (2004)  expressed  the  mass  balance  equation  of  an  aqueous-phase  constituent  in  the 
form: 


5C,  d 

-V. — ^  +  — 

dx:  dx 


a 


dC, 

dx 

UXj  J 


Q  ^sink,l  '  ^source, i 


nNAPL  _  D  dC, 
.1  '  ^source,! 

ot 


where  C/  is  the  aqueous  phase  contaminant  concentration  [ME'  ];  x/  is  distance  [L];  t  is  time  [T]; 
Vi  is  the  average  linear  groundwater  velocity  [L  T‘'];  Dy  is  the  tensor  for  the  hydrodynamic 
dispersion  coefficient  [E  T'  ];  Q  is  the  contaminant  point  source  concentration  [M  IS  ];  qs  is 
the  volumetric  flux  of  water  per  unit  volume  of  aquifer  [T'^];  6*  is  the  effective  porosity  [E'’];  Ri  is 
the  contaminant  retardation  factor  [L°];  is  a  biodegradation  mass  loss  term  dependent  on 

the  mode  of  respiration  [M  E'^  T'^];  R^Zrcei  is  a  source  term  for  the  biogenic  mass  production  [M 
E'^  T''];  and  RZt^Zi  is  a  source  term  due  to  NAPE  dissolution  [M  E'^  T'^]. 


(1) 


The  NAPE  source  term  for  petroleum  hydrocarbons  and  chlorinated  compounds  can  be 
represented  using  first-order  kinetics  for  dissolution  from  a  multiple-component  mixture 

Riss = /‘""‘-(cr  -c,)= -c,)  (2) 

where  is  the  NAPE  dissolution  rate  constant  [T‘'];  is  the  equilibrium  contaminant 

concentration  [M  E'^],  represented  as  the  product  of  the  solubility  (C™^)  of  the  constituent  in 
water  and  the  mole  fraction,  yj,  of  NAPE  constituent  /.  The  mole  fraction  is  a  function  of  the  mass 
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fractions  of  the  NAPL  constituents  and  varies  with  time.  The  mass  balance  equation  for  each 
constituent  in  the  contaminant  source  term  is  expressed  as 

r^NAPL  a 

_ _ ^  nNAPL 

T,  ^source, I 

dt  p, 

where  is  the  source  contaminant  concentration  [M  M‘^]  pb  is  the  bulk  density  of  the 

porous  medium  [Msoiid  L'^] 

For  the  multicomponent  NAPL  problems,  solutions  to  quantify  time  of  remediation  (TOR)  using 
Equations  (1)  through  (3)  are  achieved  using  numerical  models.  Sequential  Electron  Acceptor 
Model,  3D  (SEAM3D)  transport  is  a  code  designed  to  simulate  the  transport  and  attenuation  of 
contaminants  subject  to  aerobic  and  anaerobic  biotransformations  (Waddill  and  Widdowson, 
2000).  SEAM3D  includes  mass  balance  equations  for  a  multicomponent  NAPE.  The  rate  of 
release  for  each  NAPL  constituent  is  a  function  of  mass  transport  and  transfer  parameters  (v/  and 
parameters  (mass,  mole  fraction,  and  geometry),  and  chemical  properties 
(molecular  weight  and  solubility)  of  the  NAPL  components. 

This  approach  requires  an  estimate  of  the  mass  of  contaminant  present  and  an  estimate  of  the  rate 
of  ongoing  natural  attenuation  processes  acting  on  the  contaminant.  The  principal  technical 
problem,  therefore,  is  to  obtain  reliable  estimates  of  these  parameters.  Clearly,  the  reliability  of 
any  remediation  time  estimates  will  be  directly  linked  to  the  reliability  of  the  parameter 
estimates.  In  addition,  determining  remediation  times  requires  the  definition  of  an  acceptable 
contaminant  mass  threshold.  This  threshold  must  be  predetermined  in  order  to  make  remediation 
time  estimates. 

Using  solutions  to  Equation  (1),  Chapelle  et  al.  (2003)  divided  the  time  of  remediation  problem 
into  three  parts:  (1)  distance  of  plume  stabilization  (DOS),  (2)  time  of  plume  stabilization 
(TOS),  and  (3)  time  of  NAPE  dissolution  (TND).  Each  of  these  issues  can  be  addressed  using 
particular  solutions  of  Equation  (3),  which  can  be  developed  according  to  specific  needs.  NAS 
was  developed  to  make  use  of  solutions  to  Equation  (3)  in  order  to  address  these  three  classes  of 
TOR  problems.  This  software  was  designed  to  aid  the  user  in  assembling  and  organizing  the  data 
needed  to  make  TOR  estimates,  to  obtain  appropriate  and  useful  solutions  of  the  TOR  equation, 
and  to  illustrate  the  various  uncertainties  inherent  in  TOR  estimates.  No  attempt  has  been  made 
to  make  NAS  applicable  to  all,  or  even  most,  TOR  problems.  Rather,  NAS  is  designed  around 
numerous  simplifications  of  hydrologic,  microbial,  and  geochemical  processes  that,  while 
convenient,  will  introduce  unacceptable  error  to  some  problems.  The  NAS  software  can  be 
downloaded  from  the  website  http://www.nas.cee.vt.edu. 

The  NAS  tool  is  designed  for  application  to  groundwater  systems  consisting  of  relatively 
homogeneous,  saturated  media  such  as  sands  and  gravels.  In  its  present  form  NAS  is  not 
intended  for  simulating  solute  transport  in  dual-domain  porous  media  such  as  fractured  rock 
aquifers  and  highly  heterogeneous  unconsolidated  aquifers.  However,  NAS  may  be  applicable  to 
fractured-rock  aquifer  systems  which  can  be  adequately  represented  by  an  equivalent  single¬ 
porosity  system.  NAS  is  designed  for  application  to  petroleum  hydrocarbon  and  chlorinated 


4 


solvent  sites.  Version  2  of  NAS  was  recently  updated  for  an  expanded  range  of  contaminant 
groups  to  provide  greater  flexibility  for  adding  additional  contaminants  and  groups. 

2,2  PROCESS  DESCRIPTION 

First,  detailed  site  information  about  hydrogeology,  redox  conditions,  and  contaminant 
concentrations  must  be  entered.  Table  1  provides  a  summary  of  the  required  site  data.  The  goal 
of  site  data  assessment  is  to  determine  site-specific,  contaminant-specific  degradation  rates  using 
an  inverse  modeling  technique.  NAS  is  primarily  designed  as  a  screening  tool  early  in  the 
remedial  investigation/feasibility  study  (RI/FS)  process  following  completion  of  a  site 
investigation  and  characterization.  If  the  data  NAS  requires  is  not  available,  then  time  of 
remediation  estimates  cannot  and  should  not  be  made.  However,  another  use  of  NAS  is  to  reveal 
site  data  deficiencies  that  can  be  addressed  during  the  RI/FS  process  and  to  develop  monitoring 
strategies. 


Table  1,  Summary  of  NAS  v2.2  Site  Data  Requirements, 


Hydrogeology 

Hydraulic  conductivity' 

Hydraulic  gradient' 

Fraction  of  organic  matter' 

Total  porosity^ 

Effective  porosity^ 

Average  saturated  thickness  impacted  by  contamination^ 

Redox^ 

Dissolved  oxygen,  ferrous  iron,  sulfate 

Optional:  Nitrate,  Mn(II),  sulfide,  methane,  dissolved  hydrogen 

Contaminant 

Chlorinated  ethenes:  tetrachloroethene  (PCE)  (optional),  trichloroethene  (TCE),  and  daughter  products 

Petroleum  hydrocarbons:  benzene,  toluene,  ethylbenzene,  xylene,  methyl  tertiary  butyl  ether  (MTBE) 
(optional),  naphthalene  (optional) 

Chlorinated  ethanes:  trichloro ethane  (TCA)  and  daughter  products 

Chlorinated  methanes:  Carbon  tetrachloride  and  daughter  products 

Requirements 

'Best  estimate,  maximum,  minimum 

^Best  estimate 

^Values  from  3  or  more  wells  along  the  solute  plume  flow  path 

Figure  1  shows  a  flowchart  describing  how  the  NAS  software  can  be  used  to  address  TOR 
questions.  After  data  entry,  NAS  estimates  site-specific  groundwater  flow  rates,  biodegradation 
rates,  and  sorption  properties.  Based  on  the  range  of  estimates,  NAS  then  produces  either 
analytical  or  numerical  solutions  of  the  TOR  equation.  As  shown  in  Figure  1,  one  option 
employs  analytical  solutions  to  determine  the  target  reduction  in  the  source  area  concentration  to 
meet  site-specific  remediation  goals.  This  approach  and  solution  addresses  plume  concentration 
issues,  such  as  the  distance  of  stabilization  for  given  source-area  contaminant  concentrations  and 
the  time  of  stabilization  if  source-area  concentrations  are  changed.  For  the  distance  of 
stabilization,  NAS  calculates  the  allowable  maximum  source-area  concentration,  based  on  a 
regulatory  maximum  concentration  level  at  a  given  point  downgradient  of  the  source.  Then,  NAS 
estimates  how  long  it  will  take  for  the  plume  to  reach  the  lower  steady-state  configuration  once 
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source-area  eoneentrations  have  been  lowered  by  engineering  methods.  Onee  both  the  distanee 
of  stabilization  and  the  time  of  stabilization  are  acceptable,  based  on  site-speeifie  regulatory 
eriteria,  MNA  can  become  an  integral  eomponent  of  site  remediation  (Figure  1). 


Figure  1,  Flowchart  Showing  How  the  NAS  Software  Can  Be  Applied  to  TOR  Prohlems, 


The  other  option  is  a  mass-based  approaeh  to  determine  the  target  reduction  in  the  souree  area 
NAPL  or  residual  mass  to  reduee  the  TOR,  based  on  site-speeifie  remediation  goals.  To  aehieve 
this  solution  (i.e.,  when  eoneentrations  fall  below  a  given  user-supplied  threshold),  NAS  uses  the 
SEAM3D  code  (Waddill  and  Widdowson,  1998)  to  solve  Equations  (1)  through  (3)  in 
eonjunetion  with  a  groundwater  flow  eode  Modular  Three-Dimensional  Einite-Differenee 
Groundwater  Elow  Model  (MODELOW).  The  solution  provided  by  NAS  is  tailored  to  estimate 
the  length  of  time  required  by  a  given  NAPL  mass  to  dissolve  and  lower  eontaminant 
eoneentrations  at  the  souree  area.  In  prineiple,  the  numerieal  solution  could  then  be  used  to 
estimate  the  distanee  and  time  of  stabilization  for  the  remaining  residual  eoneentration.  Sinee  this 
would  signifieantly  lengthen  the  amount  of  time  required  to  eomplete  a  simulation,  numerieal 
simulation  is  not  presently  praetieal  for  distanee  and  time  of  stabilization  problems.  Rather,  once 
the  target  source  zone  eoneentration  is  determined  by  an  analytical  solution  and  the  numerical 
solution  is  completed  (time  required  to  reaeh  this  target  ealeulation),  the  analytieal  solution  for 
the  time  of  stabilization  ean  be  implemented  (Eigure  1).  Thus,  the  estimated  total  TOR  to  reaeh 
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compliance  at  a  distance  downgradient  of  the  souree  (e.g.,  plume  toe)  is  the  sum  of  the  two 
solutions. 

2.3  PREVIOUS  TESTING  OF  THE  TECHNOLOGY 

Prior  applieation  of  NAS  has  been  doeumented  in  Chapelle  et  al.  (2003)  using  limited  data  sets 
from  two  sites:  NSB  Kings  Bay,  Georgia  (ehlorinated  ethenes)  and  Marine  Corps  Air  Station 
(MCAS)  Beaufort,  South  Carolina  (petroleum  hydroearbons).  Although  these  two  sites  have 
different  eontaminant  plumes,  they  share  two  key  eharaeteristies:  (1)  high-quality,  long-term 
groundwater  monitoring  data  and  (2)  long-term  deeline  in  eontaminant  concentrations  in  souree- 
area  monitoring  wells  following  souree  remediation.  Mendez  et  al.  (2004)  presented  the 
applieation  of  NAS  at  the  MCAS  Beaufort,  South  Carolina  to  simulate  the  depletion  of  benzene 
eoneentrations  at  a  source-area  monitoring  well,  based  on  estimates  of  ground  water  veloeity  and 
light  non-aqueous  phase  liquid  (LNAPL)  mass,  eomposition,  and  dimensions. 

2.4  ADVANTAGES  AND  LIMITATIONS  OF  THE  TECHNOLOGY 

NAS  offers  several  advantages  relative  to  eomprehensive  groundwater  modeling  of  sites.  Even 
when  NAS  estimates  TOR  using  the  SEAM3D  NAPE  dissolution  option,  the  user  is  not  required 
to  speeify  numerieal  parameters  (e.g.,  grid  spacing)  or  any  spatial  input  parameters.  Beeause 
NAS  ineludes  a  simple  self-calibrating  analytical  model,  the  amount  of  time  and  effort  required 
is  mueh  less  than  for  a  site  model  for  groundwater  flow  and  solute  transport.  However,  at  the 
sites  with  eomplex  hydrogeology  and  patterns  of  groundwater  flow,  eomprehensive  groundwater 
modeling  offers  greater  capabilities  relative  to  NAS.  While  hydrogeology  and  flow  patterns 
obviously  play  a  large  role,  an  aceurate  and  complete  eharacterization  of  the  souree  area  is 
essential  to  effeetive  remediation  of  the  souree  zone  (National  Researeh  Couneil,  2004). 
Eikewise,  any  mathematieal  model  may  not  be  aecurate  if  the  souree  term  is  estimated  based  on 
limited  information.  Modeling  tools,  ineluding  NAS,  are  useful  in  developing  and  refining  site 
eoneeptual  models  of  the  aqueous  plume  and  source  zones  and  in  exposing  site  eharaeterization 
defieieneies. 
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3.0  DEMONSTRATION  DESIGN 


3.1  PERFORMANCE  OBJECTIVES 

The  objective  of  the  demonstration  is  to  assess  the  performance  of  NAS  as  a  computational  tool 
for  estimating  remediation  time  frames  following  source  zone  remediation.  To  achieve  this 
objective,  NAS  was  implemented  to  simulate  contaminant  concentrations  at  chlorinated  solvents 
sites  that  represent  a  range  of  conditions.  NAS  is  designed  to  calibrate  an  analytical  solution  to 
steady-state  plume  contaminant  concentration  data  first,  then  implement  analytical  and  numerical 
solutions  to  simulate  the  transient  response  of  plume  concentrations  due  to  source-zone 
treatment,  depletion,  or  some  combination  of  the  above. 

Several  metrics  were  used  to  evaluate  TOR  estimates:  (1)  inflection  points  of  contaminant 
concentration  profiles,  (2)  time  of  breakthrough,  (3)  slopes  of  contaminant  concentration 
profiles,  and  (4)  contaminant  concentration  profiles  for  different  NAPL  components  over  time. 
The  first  three  metrics  (inflection  points,  time  of  breakthrough,  and  slopes  of  contaminant 
concentration  profiles)  are  demonstrated  in  Figure  2,  which  depicts  total  chlorinated  ethene 
concentration  changes  at  a  location  downgradient  of  the  source  following  treatment  for  one  site 
application  (NSB  Kings  Bay,  Georgia). 


Al,  A2  =  inflection  points 
B  =  Located  at  time  of  breakthrough 
C  =  Slope  of  the  concentration  decline. 


Figure  2,  Observed  and  NAS-Predicted  Total  Chlorinated  Ethene  Concentrations 

(using  tracer-derived  ground  water  velocity  and  the  range  of  retardation  factors,  including  the 

case  (R  =  2.28)  where  error  is  minimized). 

Point  Al  represents  the  time  at  which  the  impact  of  source  remediation  is  first  observed  at  a 
downgradient  monitoring  well.  This  point  of  inflection  is  identified  at  the  time  where  a 
consistent  decline  of  contaminant  concentration  is  observed.  The  second  inflection  point,  A2 
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corresponds  to  the  time  required  to  reach  a  new  steady-state  eoneentration  at  a  specifie  loeation 
within  the  solute  plume.  Point  B  eorresponds  to  the  time  of  breakthrough  where  50%  of  the  net 
deerease  in  eoneentration  is  observed.  The  slope,  C,  is  the  rate  of  deeline  in  the  eoneentration 
between  points  A1  and  A2.  Performanee  objeetives  based  on  these  metries  are  summarized  in 
Table  2. 


Table  2.  Performance  Objectives, 


Type  of 
Performance 
Objective 

Primary  Performance 
Criteria 

Expected  Performance 
(Metric) 

Actual  Performance 
Objective  Met? 

Quantitative 

Infleetion  points  of 
eoneentration-time  data 

The  predieted  infleetion  point 
(A1  and  A2,  Figure  2)  of  the 
eoneentration-time  eurve  should 
eoineide  with  the  observed 
infleetion  point  within  one  year. 

Yes'’^’^’**;  Ineonelusive"' 

Quantitative 

Time  of  breakthrough 

The  predieted  time  of 
breakthrough  (B,  Figure  2)  of 
the  eoneentration-time  eurve 
should  eoineide  with  the 
observed  breakthrough  within  2 
years. 

Yes^’*;  Ineonelusive'  ^'’’ 

Qualitative 

Slope  of  eoneentration-time 
eurve  following  first 
infleetion  point 

The  predieted  slope  of  the 
infleetion  point  of  the 
eoneentration-time  eurve  (slope 

C,  Figure  2)  should  be  similar  to 
the  observed  slope. 

Yes"*’^’’’*;  Ineonelusive' 

Qualitative 

Predieted/Observed 
eontaminant-eoneentration 
profiles  for  different 
eomponents  of  NAPL 

NAS  will  be  assessed  for  how 
aeeurately  it  prediets  that  the 
more  soluble  NAPL  eomponents 
are  removed  sooner  than  the  less 
soluble  eomponents. 

Yes^’^;  Ineonelusive^ 

'Seneca  Army  Depot,  New  York  ^NSB  Kings  Bay,  Georgia 

^USGS  Study  Site,  New  York  "Naval  Air  Station  Cecil  Field,  Florida 

^Naval  Air  Engineering  Station  (NAES)  Lakehurst,  New  ’Naval  Air  Station  Pensacola,  Florida 

Jersey  *USGS  Study  Site,  Arkansas 

“Hill  AFB,  Utah 

3.2  SELECTION  OF  TEST  SITES 

As  a  group,  the  eriteria  and  requirements  for  site  seleetion  were  (1)  differing  hydrogeologie 
eonditions  and  (2)  various  methods  of  source  remediation.  Individually,  the  requirements  for 
sites  were  (1)  high-quality  data  sets  and  (2)  long-term  deeline  in  eontaminant  groundwater 
eoneentrations  following  souree  remediation.  Although  NAS  is  ideally  suited  for  sites  that  have 
wells  near  the  source  and  along  the  eenter  line  of  the  plume  eoupled  with  relatively  simple 
hydrogeologie  conditions,  it  was  viewed  as  useful  to  seleet  some  sites  with  non-ideal  eonditions. 
The  notion  was  not  to  set  up  the  demonstration  for  failure  but  to  test  the  limits  of  the  software 
and  assess  how  robust  the  solutions  are  in  relation  to  the  performanee  objeetives.  For  example,  at 
some  sites  the  monitoring  infrastructure  and  sample  data  resulting  from  eharacterization  efforts 
were  not  specifically  aimed  at  demonstrating  MNA  lines  of  evidenee  or  laeked  the  detail  for  a 
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comprehensive  modeling  effort.  The  eoneeptual  models  at  some  sites  were  not  eonsistent  with  a 
simple  homogeneous,  uneonsolidated  aquifer,  as  the  site  data  was  indieative  of  more  eomplex 
hydrogeologie  eonditions.  Therefore,  some  sites  seleeted  for  the  demonstration  exhibited  a  few 
non-ideal  charaeteristics. 

Desired  site  eharaeteristies  ineluded; 

•  Five  years  of  annual  monitoring  data  that  includes  the  geoehemieal  and 
eontaminant  ehemistry,  and  water  levels 

•  Appropriate  geoehemieal  data  that  allows  for  assessment  of  the  plume  status  in  an 
eleetron  aeeeptor  distribution  framework 

•  Groundwater  veloeities  exeeeding  5  ft/year 

•  Darey  flow  (i.e.,  not  eonduit  flow  sueh  as  Karst) 

•  Unidireetional  groundwater  flow  fields  and  boundary  eonditions  that  ean  be 
reasonably  simulated  by  the  analytieal  eode. 

Data  sets  from  the  seleeted  sites  demonstrated  either  natural  depletion  of  souree  zones  over  time 
or  a  deelining  trend  in  eontaminant  eoneentration  at  downgradient  wells  following  a  souree-area 
ERA.  Allowanee  was  made  for  noise  in  the  data  and/or  variation  indieative  of  sampling  or  other 
issues  that  eannot  be  sufficiently  ineorporated  into  the  solution.  Candidate  sites  were  sereened 
based  on  the  availability  of  redox  indieator  data  eolleeted  in  aeeordanee  with  EPA/DoD 
protoeols  for  MNA  (EPA,  1998).  Minimally,  dissolved  oxygen,  sulfate,  and  ferrous  iron  data 
were  neeessary.  Groundwater  velocities  and  travel  times  had  to  be  sufficient  to  measure  changes 
over  time  relative  to  the  period  of  time  for  which  monitoring  data  has  been  eolleeted.  Beeause 
the  software  is  not  designed  to  handle  eonduit  flow,  sites  with  hydrogeologie  eonditions  that  ean 
best  be  deseribed  as  eomplex  fraetured  bedroek  and/or  Karst  with  solution  ehannels  were  not 
eonsidered.  Einally,  local  groundwater  flow  pattern  along  the  length  of  the  contaminant  plume 
eould  not  be  signifieantly  impacted  by  pumping  wells  or  natural  sources  and  sinks. 

3.3  TEST  SITE,  FACILITY  HISTORY,  AND  CHARACTERISTICS 

The  eight  sites  seleeted  for  demonstrations  (Table  3)  met  the  global  eriteria  as  refleeted  in  the 
range  of  loeations,  hydrogeologie  settings,  and  remediation  options.  The  histories  of  the 
individual  sites  are  described  in  detail  in  the  Final  Report.  In  general,  the  sites  can  be 
eharaeterized  as  former  or  aetive  landfills  or  unlined  waste  disposal  faeilities.  The  hydrogeology 
at  both  sites  loeated  in  New  York  was  eharaeterized  as  non-ideal  in  which  the  flow  systems  were 
dominated  by  fractured  flow.  However,  previous  site  investigations  have  provided  reasonable 
justification  for  incorporating  Darcian  flow  in  the  site  conceptual  models.  Uniform  groundwater 
flow  field  was  a  eommon  eharaeteristic  at  all  sites  with  the  exeeption  of  two  oases — the  pump- 
and-treat  site  and  Hill  AFB.  At  the  latter,  the  presenoe  of  the  barrier  wall  impeded  flow,  resulting 
in  a  nonuniform  horizontal  hydraulio  gradient.  Furthermore,  the  eontaminant  plume  is  present 
within  two  formations,  whioh  results  in  a  ohange  in  hydrogeologie  parameters  along  the 
fiowpath. 


11 


Each  site  also  met  the  two  requirements  for  data  quality  and  deelining  eoneentration  time  trends. 
Chlorinated  ethenes  were  a  eommon  eontaminant  to  the  sites,  whieh  ineluded,  but  were  not 
limited  to,  PCE,  TCE,  diehloroethene  (DCE)  and  vinyl  ehloride  (VC).  Chlorinated  benzenes 
were  present  at  both  sites  loeated  in  Elorida. 

Table  3,  Summary  of  Demonstration  Sites, 


Facility 

Site 

Hydrogeology 

Source  Zone  ERA 

Seneca  Army  Depot,  NY 

Ash  Landfill 

Glacial  till  and  weathered 
shale 

Excavation 

USGS  Study  Site,  NY 

Textron 

Fractured  dolomite 

Pump-and-treat 

NAES  Lakehurst,  NJ 

Sites  I  &  J 

Fine  to  coarse  quartz  sand 

Natural  depletion 

Hill  AFB,  UT 

Operating  Unit  2 
(OU2) 

Fluvial  fine  sand  and  silt  and 
interbedded  clay  layers 

Containment  (barrier 
wall) 

NSB  Kings  Bay,  GA 

Site  11 

Medium  sand  (marine) 

Chemical  oxidation 
(Fenten’s) 

Naval  Air  Station  Cecil  Field, 
FL 

Site  3 

Fine  to  medium  quartz  sand 
&  interbedded  clay  and  sandy 
clay  layers 

Air  sparging 

Naval  Air  Station  Pensacola, 
FL 

WWTP 

Fine  to  medium  quartz  sand 
(marine  and  fluvial) 

Oxygen  release 
compound 

USGS  Study  Site,  AK 

Alaska  Department 
of  Transportation 
(DOT) 

Highly  permeable  alluvium 

Excavation 

3.4  PHYSICAL  SETUP  AND  OPERATION 

Physieal  setup  and  operation  of  NAS  is  deseribed  in  Seetion  2.2,  Proeess  Deseription.  Eor  the 
site  demonstrations,  NAS  was  implemented  in  a  predietive  fashion  using  the  site-specifie  ranges 
in  hydrogeologie,  sorption,  and  geoehemieal  (redox)  parameters.  A  eomparison  of  the  observed 
data  (long-term  performanee  monitoring)  to  the  range  in  simulated  eoneentration  versus  time 
proeeeded  without  ealibration  or  adjustment  of  parameters. 

3.5  SAMPLING/MONITORING  PROCEDURES 

Not  applieable  to  the  demonstration. 

3.6  ANALYTICAL  PROCEDURES 
Not  applieable  to  the  demonstration. 
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4.0  PERFORMANCE  ASSESSMENT 


4.1  PERFORMANCE  CRITERIA 

All  eight  site  demonstrations  are  described  in  detail  in  the  Final  Report  to  ESTCP  (2006).  Below 
is  a  summary  of  the  simulation  results  in  comparison  to  the  observed  data. 

Time  of  stabilization  estimates  were  compared  to  observed  data  at  monitoring  wells  located 
downgradient  of  the  source  at  distances  ranging  from  140  to  716  ft  from  the  source  without 
adjustment  of  input  parameters  (i.e.,  direct  comparison  of  observed  versus  simulated  with  no 
calibration).  Applications  of  NAS  at  both  at  the  Seneca  Army  Depot  Activity  (SEDA)  Ash 
Landfill  and  Hill  AEB  OU2  (Eigures  3  and  4,  respectively)  demonstrate  that  fluctuations  over 
time  in  the  magnitude  and  direction  of  the  hydraulic  gradient  will  impact  the  position  of  the 
plume  centerline  and  will  contribute  to  concentration  variability  relative  to  the  ideal  solution. 

The  response  in  TCE  concentration  at  two  monitoring  wells  at  the  SEDA  Ash  Landfill  is  shown 
in  Eigure  3  following  source  excavation.  Data  gaps  in  the  mid-1990s  and  the  time  fluctuation  in 
concentrations  hinder  the  accurate  estimation  of  the  observed  breakthrough  time.  The 
breakthrough  and  inflection  points  are  difficult  to  judge  from  the  data  but,  in  general,  a 
reasonable  fit  was  obtained  at  both  wells.  Eor  both  wells,  the  TOS  estimates  are  consistent  with 
the  observed  data.  The  likely  cause  of  data  variability  is  the  reduction  in  the  source  volume 
through  excavation  combined  with  fluctuations  in  the  direction  of  the  horizontal  hydraulic 
gradient  with  time.  The  non-ideal  response  can  also  be  attributed  in  part  to  variability  in  the 
groundwater  flow  path  through  the  weathered  shale  and  till. 


TCE  .It  Well  PT12A  .it  200.00  ft  TCE  .it  Well  PT21A  .it  400.00  ft 


Figure  3,  Observed  and  NAS-Simulated  TCE  Concentration 
Versus  Time  at  Downgradient  Wells 

(based  on  the  observed  source  concentration  reduction  at  the  SEDA  Ash  Landfill). 


Analytical  solutions  used  in  NAS  are  based  on  the  assumption  that  the  hydrostratigraphic  unit  is 
homogeneous  so  that  contaminant  velocities  do  not  vary  along  the  entire  length  of  the  plume.  At 
this  site  (OU2)  the  plume  is  located  in  two  different  aquifer  systems.  In  this  application. 
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parameters  are  ealeulated  based  on  the  parameters  of  the  two  separate  formations.  Figure  4 
shows  that  simulations  based  on  the  range  of  parameters  derived  from  the  formation  closest  to 
the  containment  wall  resulted  in  the  closest  match  with  the  observed  data  at  both  wells.  Although 
NAS  was  not  designed  to  simulate  the  effects  of  a  containment  wall  on  concentrations,  these 
results  suggest  that  NAS  is  useful  in  simulating  the  plume  changes  over  time. 


TCE  at  Well  U2-021R  at  407.00  ft 


TCE  at  Well  U2-043  at  716.50  ft 


Time  [yr] 


Figure  4,  TCE  Concentration  Versus  Time  at  Wells  Downgradient  of  the  Slurry  Wall  at 

Hill  AFB  OU2 

(Contaminant  velocity  range  is  based  on  hydrogeological  parameters  derived  from  the  near¬ 
source  hydrostratigraphic  unit  and  a  site-averaged  hydraulic  gradient.) 

A  third  application  (Figure  2)  utilized  a  tracer-based  value  of  the  groundwater  velocity.  A  pulse 
of  sulfate  was  observed  in  several  observation  wells  resulting  from  the  use  of  Fenton’s  reagent  in 
the  source  zone  and  subsequent  transport  with  the  natural  gradient.  The  results  using  a  tracer- 
based  velocity  estimate  showed  improved  agreement  between  the  observed  and  calculated 
concentration  trend  when  compared  to  the  initial  estimates. 

The  NAS  Cecil  Field,  Site  3  demonstration  illustrates  the  fourth  performance  objective  involving 
the  multi-component  NAPL  dissolution  concept  and  estimates  for  TND  in  the  source  area.  Figure 
5  depicts  concentration  versus  time  curves  for  six  NAPL  components  relative  to  historical  data, 
using  a  range  of  groundwater  velocities.  The  analysis  suggests  that  the  volatile  organic 
compound  (VOC)  concentration  trends  during  the  early  stages  of  NAPL  dissolution  appear  to  be 
most  sensitive  to  estimates  of  mass.  During  these  periods  the  more  soluble  compounds  elute 
early  and  a  more  pronounced  peak  is  observed  in  the  plume  concentrations.  Because  the 
observed  data  more  closely  follows  the  average  estimate  (5000  kg)  of  mass  for  both  VOC  and 
semi-volatile  organic  compounds  (SVOC),  this  value  was  used  to  estimate  the  TND.  The 
potential  range  of  transport  velocities  predicts  a  range  of  values  of  TND  that  encompasses  the 
observed  data.  The  results  demonstrate  that  the  SEAM3D  NAPL  package  is  robust  and  capable 
of  simulating  the  differing  concentration  changes  over  time,  depending  on  the  solubility  and 
other  physical-chemical  properties  of  the  NAPL  components.  This  also  illustrates  the  combined 
effect  of  source  remediation  with  MNA. 
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Concentration  [|jgjL]  Concer^ration  (pg.!]  Concerrtration  (pg.!] 


12DCB 


TCE 


Figure  5,  Observed  Versus  NAS-Simulated  Concentrations  at  a  Source-Zone  Monitoring 
Well  for  Six  NAPL  Compounds  Using  a  Range  of  Velocity  Estimates, 

4.2  PERFORMANCE  CONFORMATION  METHODS 

Performance  criteria  and  actual  performance  of  NAS  are  summarized  in  Table  4. 
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Table  4,  Performance  and  Performance  Confirmation  Methods. 


Performance  Criteria 

Expected 

Performance  Metric 
(pre-demonstration) 

Performance 
Confirmation  Method 

Actual  (post¬ 
demonstration) 

Primary  Criteria  (Performance  Objectives) 

(Quantitative) 

Inflection  point  of 
concentration-time 

curve 

The  predicted  inflection 
point  (point  At,  Fig.  2)  of 
the  concentration-time 
curve  should  coincide 
with  the  observed 
inflection  point  within 
one  year. 

Concentration-time 
curves  at  least  one  mid¬ 
plume  well,  and 
preferably  several,  over 
time 

In  general,  better 
performance  noted  at 
wells  closer  to  source 
relative  to  further 
downgradient  wells 

Time  of  breakthrough 

The  predicted  time  of 
breakthrough  (point  B, 

Fig.  2)  of  the 
concentration-time  curve 
should  coincide  with  the 
observed  breakthrough 
within  two  years. 

Concentration-time 
curves  at  least  one  mid¬ 
plume  well,  and 
preferably  several,  over 
time 

NAS  captured  velocity 
uncertainty  in  time  of 
breakthrough. 

Performance  was  superior 
at  sites  where  flow 
velocity  is  constrained 
through  independent  data. 

Primary  Criteria  (Performance  Objectives) 

(Qualitative) 

Slope  of  couceutratiou- 
time  curve  followiug 
iuflectiou  poiut 

The  predicted  slope  of  the 
inflection  point  of  the 
concentration-time  curve 
(slope  C,  Figure  2)  should 
be  similar  to  the  observed 
slope. 

Concentration-time 
curves  at  least  one  mid¬ 
plume  well,  and 
preferably  several,  over 
time 

Greater  data  variability 
due  to  fluctuations  in 
hydraulic  gradient  or 
other  factors 
downgradient  wells 

Predicted/observed 
coutamiuaut- 
couceutratiou  profiles 
for  differeut 
compoueuts  of  NAPL 

More  soluble  components 
of  NAPL  should  be 
removed  sooner  than  less 
soluble  components 
relative  to  overall  NAPL 
composition. 

Concentration-time 
curves  at  near-source 
well,  and  preferably 
several,  over  time 

Adequate  performance  for 
the  data  time  windows 
(10-years)  considered. 

4,3  DATA  ASSESSMENT 

Section  4.1,  Performance  Criteria,  contains  a  detailed  assessment  of  the  NAS  results.  The 
following  summarizes  the  findings  with  respect  to  accuracy,  versatility,  reliability,  and 
applicability. 


•  Overall,  the  performance  of  NAS  was  satisfactory  in  matching  the  inflection 
points  and  time  of  breakthrough  for  the  five  sites  where  TOS  was  evaluated. 

•  Using  the  field-measured  source  concentration  following  remediation  and  the  site- 
specific  estimated  contaminant  velocities,  the  best  match  between  the  observed 
and  simulated  concentrations  was  achieved  at  the  monitoring  wells  closest  to  the 
source. 
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•  Less  accurate  results  were  observed  in  some  cases  at  wells  further  downgradient 
at  monitoring  wells  due  to  time  fluctuations  in  the  direction  and  magnitude  of  the 
groundwater  velocity. 

•  NAS  was  capable  of  simulating  the  dissolution  and  dissipation  of  individual 
NAPL  components  at  the  three  sites  where  TND  simulations  were  performed. 

•  The  hydrogeology  of  the  sites  did  not  appear  to  be  a  factor  in  the  performance  of 
NAS  at  any  of  the  eight  applications. 

•  The  mode  of  source  remediation  was  not  a  factor  in  the  performance  of  NAS. 

•  No  significant  reliability  problems  were  encountered  in  the  implementation  of 
NAS  version  2.2.0  at  any  of  the  eight  sites. 

•  The  results  of  this  demonstration  strongly  suggest  that,  because  NAS  is  based  on 
sound  science,  it  can  serve  as  an  effective  tool  for  decision  making,  data  analysis, 
and  cost  optimization  at  a  wide  range  of  contaminated  sites  and  is  not  limited  to  a 
small  subset  of  “simple  sites”  because  of  its  simplicity. 

4,4  TECHNOLOGY  COMPARISON 

The  demonstration  of  NAS  at  NAES  Lakehurst,  New  Jersey,  Sites  I  &  J  provided  a  comparison 
of  technologies;  in  this  case  with  a  comprehensive  numerical  solute  transport  model  constructed 
using  the  code  RT3D  and  also  using  a  regression  analysis  with  extrapolation  forward  to 
determine  TOR.  NAS  was  employed  to  predict  the  decline  in  aqueous  concentration  resulting 
from  the  dissolution  of  the  NAPL  source  at  a  nearby  downgradient  monitoring  well.  Because  the 
composition  and  location  of  the  source  could  not  be  well-defined  during  the  monitoring  program, 
two  hypothetical  source  compositions  were  modeled  using  NAS:  (1)  single-composite  NAPL 
source  (total  chlorinated  ethenes)  and  (2)  multiconstituent  NAPL.  The  comparison  is  provided 
for  the  first  case  only. 

Using  NAS  and  calibrating  to  the  first  5  years  of  data  only,  the  single  compound  assumption 
yielded  a  best  fit  to  the  observed  data  using  a  source  mass  of  8,000  kg  and  a  mass  fraction  of 
0.05  (Figure  6).  Unlike  SEAM3D,  which  employs  a  mass  balance  approach  to  the  source  term, 
depleting  sources  in  RT3D  must  be  manipulated  by  the  model  users  through  the  implementation 
of  the  model  input  (e.g.,  time-varying  boundary  condition).  Source  mass  depletion  must  be 
calculated  outside  of  RT3D  (typically,  in  a  spreadsheet)  so  that  TOR  is  based  on  an  assumed 
starting  mass  and  model-simulated  mass  flux  for  the  rate  of  depletion.  The  model  user  must 
make  assumptions  about  the  relationship  between  source  mass  and  contaminant  concentrations  in 
groundwater  (e.g.,  directly  correlated)  for  implementation  of  the  site  model.  The  failure  to 
capture  the  time  trend  using  this  approach  is  evident  in  the  RT3D  results  (Figure  6).  The  simple 
data  extrapolation  method  was  employed  using  the  complete  data  set.  The  inherent  assumption 
with  this  approach  is  that  source  decay  (i.e.,  groundwater  concentration  decline)  is  first  order  and 
will  continue  to  decay  at  the  same  rate.  A  comparison  with  the  source  mass-balance  approach 
suggests  that  a  first-order  solution  may  oversimplify  the  problem,  resulting  in  overly  optimistic 
estimates  of  TOR. 
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Total  Chlorinated  Ethene  Concentration  (|ag/L) 


Figure  6,  Observed  Concentrations  of  Total  Chlorinated  Ethenes  for  Single-Component 

Source  Models  Using  NAS  and  RT3D 

(A  first-order  regression  model  is  also  shown.) 
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5.0  COST  ASSESSMENT 


5.1  COST  REPORTING 

The  NAS  software  is  used  to  estimate  the  time  in  whieh  remedial  objeetives  will  be  met.  Once 
this  time  frame  is  established,  the  cost  to  monitor  the  site  can  be  estimated.  Typical  costs 
associated  with  a  long-term  performance  monitoring  program  include  the  cost  of  sampling  and 
analysis  of  wells  and  an  annual  report  compiling  the  progress  of  remediation  for  the  site.  Often, 
long-term  monitoring  programs  incorporate  more  frequent  sampling  at  the  beginning  of  the 
monitoring  program  and  a  less  frequent  sampling  as  trends  become  established  and  seasonal 
impacts  are  understood.  At  sites  where  source  remediation  is  combined  with  MNA,  the  cost  of 
implementing  source  ERA  is  included.  For  this  demonstration,  the  estimated  cost  was 
determined  to  sample  the  number  of  wells  at  the  NAS  Cecil  Field  site  until  the  estimated  time  to 
reach  remedial  objectives  was  met.  The  basis  for  the  cost  estimate  is  summarized  in  Table  5. 

NAS  may  be  used  as  part  of  an  RI/FS  and  in  a  postaudit  analysis  of  a  site  where  MNA  with  or 
without  source  remediation  is  the  selected  remedy.  Under  both  of  these  scenarios  the  time  and 
cost  to  build  a  comprehensive  flow  and  transport  numerical  model  can  require  significant  effort 
including  the  additional  cost  associated  with  required  field  data.  For  this  demonstration  the 
estimated  cost  to  implement  NAS  was  determined  and  compared  to  the  estimated  cost  associated 
with  implementing  a  comprehensive  flow  and  transport  model  to  evaluate  MNA.  The  analysis 
compares  the  cost  of  using  the  two  technologies  early  in  the  RI/FS  with  the  same  data 
availability. 


Table  5,  Basis  for  Cost  Estimates, 


Cost  Factor 

Variables 

Number  of  wells  in  sampling  program 

Optimized  sampling  program? 

Frequency  of  sampling 

Optimized  sampling  program? 

Sample  collection 

Geographic  labor  rates  and  well  characteristics 

Field  and  laboratory  analysis 

Fixed  base  laboratory  and  field  portable  kits 

Time  to  meet  remedial  objectives 

Estimated  within  a  range  by  NAS 

Source  zone  removal 

Yes/no? 

Annual  reporting 

Geographic  labor  rates 

5.2  COST  ANALYSIS 
5.2.1  Life-Cycle  Cost  Analysis 

An  estimate  of  the  life-cycle  costs  of  monitoring  associated  with  TND  for  the  current  MNA 
remedy  at  NAS  Cecil  Field,  Site  3  was  developed  in  conjunction  with  NAS  TOR  estimates.  The 
present  cost  of  semi-annual  monitoring  for  VOC,  SVOC,  and  geochemical  parameters  at  11 
monitoring  wells  is  shown  in  Table  6.  The  cost  of  preparing  an  annual  report  is  $100,000.  The 
net  present  value  of  long  term  monitoring  at  Site  3  is  estimated  in  Table  6  using  both  the 
estimated  TND  for  TCE  (43  years)  and  naphthalene  (>69  years)  to  achieve  NFA  criteria.  The  net 
present  value  for  monitoring  over  this  period  of  time  is  at  least  $2,333,029.  The  accuracy  of 
these  cost  estimates  is  largely  dependent  on  the  TOR  and  monitoring  requirements.  As  a  result  of 
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the  high  TOR  estimates  and  assoeiated  life  eyele  costs  as  well  as  the  loss  of  efficiency  associated 
with  MNA  of  the  chloroethenes  in  the  source  area  resulting  from  air  sparging,  additional 
optimization  at  Site  3  is  recommended. 


Table  6.  Estimated  Cost  for  MNA  at  NAS  Cecil  Field,  Site  3. 


Description 

Annual  Costs 

Performance  monitoring,  reporting,  and  project 
management 

$100,000 

Total  (TCE)‘ 

$2,037,079^ 

Total  (Naphthalene)^ 

$2,333,029^ 

'TND  =  43  years 
^TND  =  >69  years 

^Assuming  a  4%  return  on  investment 


5.2.2  Implementation  Cost  Analysis 

An  estimate  of  the  cost  to  implement  a  model  in  conjunction  with  the  scope  of  work  associated 
with  the  RI/FS  is  highly  site-specific.  Comprehensive  models  for  groundwater  flow  and 
contaminant  transport  are  often  used  to  evaluate  a  range  of  remedial  strategies.  Direct  cost 
comparison  to  a  previously  implemented  site  (e.g.,  NAES  Lakehurst)  is  problematic  because 
work  associated  with  implementation  of  models  is  typically  buried  inside  the  total  cost  of 
conducting  a  feasibility  study  at  a  site. 

Table  7  provides  a  comparison  of  implementation  costs  for  the  two  methods,  quantifying  TOR 
using  NAS  (v2)  and  using  a  comprehensive  model.  Although  NAS  is  designed  to  be  used  by 
personnel  who  may  not  have  adequate  education  and/or  training  to  implement  a  comprehensive 
model,  the  same  hourly  rate  ($  150/hr)  is  used  for  this  comparison.  The  development  of  a 
modeling  plan  entails  evaluation  of  site  data  (hydrogeologic,  contaminant,  and  redox  indicator 
data),  description  of  a  conceptual  model,  and  detailed  plans  for  construction  of  the  numerical 
model  and  parameter  estimation.  The  latter  task  is  not  required  for  NAS  but  constitutes  a 
significant  labor  cost  for  a  comprehensive  model.  As  reflected  in  the  costs  shown  in  Table  7,  less 
documentation  is  required  for  NAS.  To  implement  a  flow  and  transport  model  at  the  same  site, 
larger  time  requirements  are  largely  due  to  the  efficiency  of  NAS  and  the  difficulty  in 
constructing  and  calibrating  a  fully  3-D  (or  even  2-D)  numerical  model.  In  addition,  the  reporting 
requirements  are  substantially  greater  for  documenting  the  comprehensive  approach. 

The  ratio  of  the  estimated  cost  of  using  NAS  compared  to  the  cost  of  using  a  comprehensive 
model  is  equal  to  the  ratio  of  the  estimated  hours  (5.6:1)  to  complete  the  analysis  ($19,200  and 
$108,000,  respectively).  Note  that  the  costs  shown  in  this  comparison  do  not  account  for 
establishment  of  a  contract,  management  of  this  contract,  or  contracting  rates  and  fees.  It  is 
reasonable  to  assume  that  the  sum  of  these  costs  will  be  larger  for  an  investigation  using  a 
comprehensive  model.  In  addition,  it  is  reasonable  to  assume  that  level  of  oversight  by  both  a 
senior  technical  expert  and  regulating  agencies  will  be  greater  in  this  case  relative  to  NAS, 
resulting  in  higher  costs  for  the  use  of  a  comprehensive  model.  Thus,  the  cost  comparison 
(Table  7)  may  represent  only  a  third  or  less  of  the  overall  cost. 
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Table  7.  Cost  Comparison  of  TOR  Analysis  Using  NAS  Versus  a  Comprehensive  Model, 


Task  Description 

NASv2 

Comprehensive  Model 

Modeling  plan 

$6,000 

$31,500 

Model  eonstmetion 

$1,200 

$12,000 

Calibration  and  sensitivity  analysis 

$2,400 

$24,000 

Tor  simulation 

$3,600 

$22,500 

Reporting 

$6,000 

$18,000 

Total 

$19,200 

$108,000 

Another  cost-saving  feature  of  NAS  is  the  number  of  monitoring  wells  required  for  TOR 
analysis.  In  the  comprehensive  modeling  approach,  data  (concentration  and  hydraulic  head)  from 
monitoring  wells  throughout  the  site  are  used  for  model  calibration.  Because  NAS  only  requires 
data  from  the  approximate  plume  centerline,  costs  associated  with  sample  collection  and  analysis 
can  be  substantially  reduced  using  NAS. 

For  purpose  of  a  comparison  of  the  total  cost  savings,  two  hypothetical  sites  are  considered;  the 
“small  site”  and  the  “large  site,”  consisting  of  20  and  50  monitoring  wells,  respectively.  The 
small  site  consists  of  a  stable  attenuated  plume  <100  m  in  length.  At  the  large  site,  the 
combination  of  the  source  mass  flux  and  natural  attenuation  capacity  results  in  a  longer,  but 
stable,  plume.  Although  the  number  of  active  monitoring  wells  varies  considerably  from  site  to 
site,  these  values  are  for  the  purpose  of  illustration  and  calculation  of  cost  savings.  In  the  case  of 
the  small  site,  it  is  assumed  that  seven  monitoring  wells  will  reasonably  represent  the  plume 
centerline.  For  the  large  site,  the  number  of  wells  required  to  define  the  centerline  is  assumed  to 
be  eleven.  These  assumptions  are  based  in  part  on  the  NAS  site  applications  described  in  this 
report.  Costs  for  laboratory  and  field  analysis  of  groundwater  constituents  used  for  this 
comparison  are  $300  per  sample  for  quantifying  contaminant  concentrations  and  $150  per 
sample  for  redox  indicator  data  (inorganic  constituents  and  dissolved  oxygen).  Analytical  costs 
are  for  a  one-time  event  and  do  not  reflect  monitoring  over  the  life  cycle  of  a  site.  Costs  to 
implement  a  modeling  approach  at  small  sites  are  given  in  Table  7.  For  the  large  site  cost 
analysis,  the  increases  in  modeling  costs  equal  10%  and  20%  for  the  small  and  large  sites, 
respectively. 

Table  8  summarizes  the  total  cost  for  each  case  using  the  two  modeling  approaches.  Analytical 
cost  components  are  based  on  the  rates  presented  above,  while  modeling  cost  components  are 
presented  in  Table  7.  The  results  show  an  overall  per  site  cost  savings  between  $94,650  and 
$126,030  can  be  realized  by  using  NAS  as  the  modeling  approach.  Although  cost  savings  at  any 
given  site  may  exceed  or  be  less  than  this  range,  the  results  suggest  that  considerable  savings  can 
be  achieved  to  calculate  TOR  using  NAS  relative  to  comprehensive  modeling  (Figure  7).  Since 
sampling  labor  was  not  accounted  for,  per  application  cost  savings  could  be  considerably  higher. 
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Table  8,  Cost  Comparison  and  Savings  for  Two  Hypothetical  Sites  of  Differing  Size, 


Site 

Conventional 

NASv2 

Cost  Savings 

Samples 

Analytical^ 

Totaf’^ 

Samples 

Analytical* 

Totaf’^ 

Analytical 

Total 

Small 

20 

$9,000 

$117,000 

7 

$3,150 

$22,350 

$5,850 

$94,650 

Large 

50 

$22,500 

$152,100 

11 

$4,950 

$26,070 

$17,550 

$126,030 

‘Analytical  sampling  and  analysis  costs  =  $300/well  for  contaminants  and  $  150/well  for  redox  indicators 
"Total  costs  =  Sum  of  analytical  and  modeling  (Table  7) 

^Cost  increases  for  modeling  large  sites  =  20%  and  10%  using  conventional  models  and  NAS,  respectively 


Figure  7,  Cost  Comparison  for  TOR  Analysis  at  a  Single  Site  Using  a  Comprehensive 
Model  and  NAS  for  Two  Cases:  Small  and  Large  Sites 

(20  and  50  monitoring  wells,  respectively). 

A  cost  comparison  illustrating  potential  cost  savings  achieved  at  multiple  sites  over  a  5-year 
period  is  shown  in  Figure  8.  The  cumulative  cost  savings  potentially  realized  by  implementation 
of  NAS  on  a  national  level  by  DoD  RPM  was  calculated  by  assuming  that  NAS  is  being  used  to 
calculate  TOR  on  an  incremental  basis  over  time.  In  this  scenario,  one  site  per  month  is  modeled 
using  NAS  in  year  1.  In  years  2  and  3,  the  number  of  sites  increases  to  two  and  four  per  month, 
respectively,  as  a  result  of  technology  transfer  efforts  focused  on  the  use  and  benefits  of  NAS.  In 
years  4  and  5,  the  number  of  sites  evaluated  using  NAS  becomes  five  per  month,  resulting  in  the 
total  of  204  sites  evaluated  using  NAS.  During  the  5 -year  period,  it  is  assumed  that  75%  of  the 
NAS  applications  take  place  at  small  sites. 

The  cumulative  5-year  cost  savings  to  DoD  by  implementing  NAS  relative  to  a  comprehensive 
model  at  contaminated  sites  for  the  purpose  of  estimating  TOR  is  $20,900,000.  Because  this 
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comparison  does  not  factor  in  a  rate  of  savings  over  time,  the  total  may  underestimate  the  true 
savings.  This  analysis  does  not  factor  in  the  additional  potential  cost  savings  to  DoD  at  these 
sites  by  implementing  an  MNA-based  strategy  with  or  without  source  zone  reduction.  This  type 
of  analysis  is  highly  site-specific,  as  illustrated  by  the  NAS  Cecil  Field  Site  3  example  described 
in  Section  5.2.1,  Life-Cycle  Cost  Analysis. 


Figure  8,  Cost  Savings  Using  NAS  Over  a  5-Year  Period 

(beginning  with  12  sites  in  the  first  year,  24  sites  in  the  second  year,  and  48  sites  in  each  of  the 
remaining  3  years.  The  ratio  of  small  sites  to  large  sites  is  4:1). 
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6.0 


IMPLEMENTATION  ISSUES 


6.1  COST  OBSERVATIONS 

NAS  is  designed  for  use  by  site  managers  and  eonsultants  who  have  an  educational  background 
in  groundwater  hydrology  and  contaminant  transport  but  who  do  not  necessarily  have  advanced 
degrees  in  the  area  of  numerical  modeling.  The  cost  to  contract  a  consultant  to  implement  NAS 
is  consistent  with  a  consultant’s  hourly  rate  to  conduct  an  investigation  using  a  comprehensive 
numerical  model  (with  an  identical  TOR-based  modeling  objective).  The  primary  difference  in 
cost  is  attributed  to  the  greater  number  of  hours  required  for  the  latter  to  develop  and  set  up  the 
site  model  (both  groundwater  flow  and  solute  transport)  to  perform  calibration  and  sensitivity 
analyses  and  to  document  the  work.  Given  a  common  starting  point  where  data  needs  are 
sufficient  for  both  approaches,  it  is  estimated  that  the  time  required  to  complete  the  analysis  with 
NAS  is  no  more  than  20%  of  the  total  number  of  hours  expected  to  complete  an  investigation 
with  a  comprehensive  flow  and  transport  model.  It  is  important  to  note  that  one  conclusion  of  an 
application  of  NAS  is  the  need  to  apply  a  comprehensive  numerical  modeling  to  a  site  to  more 
adequately  address  site  remediation  options  and  TOR  questions.  However,  the  initial  work  with 
NAS,  including  data  assembly  and  model  results,  will  likely  benefit  the  additional  work  and  lead 
to  cost  savings  through  efficiency  of  effort. 

6.2  PERFORMANCE  OBSERVATIONS 

NAS  met  all  performance  objectives  and  criteria.  In  cases  where  significant  model  assumptions 
were  not  violated,  NAS  was  able  to  provide  accurate  estimates  in  the  range  of  TOR.  For  the 
particular  case  of  TOS  calculations,  NAS  was  applied  in  a  predictive  mode  without  calibration. 
Thus,  the  performance  assessment  of  NAS  to  match  observed  data  was  not  based  on  the  results 
of  a  data  fitting  exercise.  In  addition  to  accuracy,  the  versatility  of  NAS  was  demonstrated  in  a 
variety  of  hydrogeological  settings.  Four  sites  were  located  in  coastal  plain  regions,  and  the  other 
four  sites  reflect  more  diverse  and  complex  hydrogeologic  settings.  For  the  latter  set  of  sites,  a 
reasonable  match  between  observed  data  and  NAS  simulation  results  of  TOS  and  TND  was 
obtained.  NAS  was  employed  at  sites  where  several  different  source  remediation  technologies 
were  applied.  NAS  proved  to  be  reliable  and  applicable  to  each  of  the  eight  sites.  Even  for  the 
case  of  a  containment  wall  (Hill  AFB)  where  the  assumption  of  uniform  flow  was  violated,  the 
NAS  application  resulted  in  an  adequate  match  with  observed  data  at  monitoring  wells  located  a 
sufficient  distance  away  from  the  wall  where  the  groundwater  flow  returned  to  a  natural 
condition. 

6.3  SCALE-UP 

Not  applicable. 

6.4  OTHER  SIGNIFICANT  OBSERVATIONS 

Performance  of  NAS  to  predict  the  equilibrium  concentration  at  a  downgradient  monitoring  well 
following  source  remediation  (TOS)  is  dependent  on  data  quality  (in  this  case,  concentration  of 
VOC  in  groundwater)  and  the  proximity  of  the  monitoring  to  the  plume  centerline.  Experience 
has  shown  that  data  assembly  and  the  development  of  a  site  conceptual  model  are  not  trivial 
efforts.  It  is  beneficial  to  complete  this  step  before  proceeding  with  an  application  of  NAS  to  a 
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site.  The  data  requirements  for  NAS  are  not  excessive  but  should  be  considered  in  data  collection 
strategies. 

6.5  LESSONS  LEARNED 

NAS  can  be  applied  to  sites  either  (1)  as  a  stand-alone  screening  tool  using  preremediation  data 
to  determine  what  contaminant  concentration  leaving  the  source  zone  is  needed  to  meet  a  site- 
specific  remedial  action  objective  (RAO)  or  (2)  following  source  remediation  using  observed 
source  zone  data.  For  this  demonstration,  NAS  was  employed  for  the  latter  case  where  the 
observed,  postremediation  concentration  at  monitoring  wells  immediately  downgradient  of  the 
source  zone  served  as  input  to  NAS. 

In  testing  the  DOS/TOS  option,  the  point  of  the  application  at  five  sites  was  to  show  what  NAS 
predicts  at  downgradient  wells  without  fitting  the  results  to  the  data  (i.e.,  it  was  not  a  curve¬ 
fitting  exercise).  The  NAS  applications  honored  the  source  concentration  following  remediation 
observed  at  the  field  sites,  and  these  inputs  to  NAS  were  not  tweaked  or  modified  to  improve  the 
fit  to  the  data.  In  either  case,  NAS  self-calibrates  a  steady-state  solution  to  the  preremediation 
plume  concentration  data  and  then  calculates  concentration  versus  time  at  downgradient 
monitoring  wells  using  the  field-measured  source  concentration  following  remediation.  The  best 
matches  between  the  observed  and  simulated  concentrations,  including  inflection  points,  were 
achieved  at  the  monitoring  wells  closest  to  the  source  (travel  distances  ranging  from  1 8  to  407  ft 
with  an  average  of  185  ft  from  source).  In  general,  less  accurate  results  were  observed  at 
monitoring  wells  located  further  downgradient  from  the  source.  At  one  site  where  excavation 
was  employed  for  source  remediation  (Seneca),  a  decrease  in  the  apparent  source  width  was 
observed  following  remediation  that  repositioned  the  plume  centerline.  NAS  accounts  only  for  a 
reduction  in  the  contaminant  concentration  leaving  the  source  zone  and  does  not  account  for 
changes  in  the  source  width  following  remediation.  The  issue  is  amplified  at  sites  where  time 
fluctuations  in  the  direction  and  magnitude  of  the  groundwater  velocity  impact  plume 
concentrations.  Time  fluctuations  in  concentration  can  be  expected  in  situations  where  the 
contaminant  mass  flux  “misses”  the  well  if  the  direction  of  groundwater  flow  changes 
seasonally.  The  most  accurate  results  were  observed  at  NSB  Kings  Bay  where  the  groundwater 
velocity  was  calculated  based  on  tracer  test  data.  These  findings  suggest  that  accurate 
determinations  of  the  groundwater  velocity  and  sorption  parameters  are  required  for  improved 
accuracy  of  TOS  following  source  zone  remediation  at  all  sites. 

The  most  accurate  results  for  concentration  reductions  following  source  remediation  (TOS)  were 
observed  at  the  NSB  Kings  Bay,  Georgia,  site  where  the  groundwater  velocity  was  calculated 
based  on  tracer  test  data.  These  results  suggest  that  because  of  uncertainty  regarding 
hydrogeology,  that  NAS-derived  TOR  results  (and  not  unlike  any  other  model  results)  should  not 
be  viewed  as  a  precise  answer  but  as  an  initial  estimate  that  will  need  further  refining  over  time. 
These  findings  suggest  that  site  resources  directed  toward  quantifying  groundwater  velocity  and 
sorption  parameters  will  reduce  uncertainly  of  TOS  estimate  following  source  zone  remediation. 

Results  of  the  TND  simulations  show  that  NAS  was  capable  of  simulating  the  dissolution  and 
dissipation  of  individual  NAPL  components  at  the  remaining  three  sites  listed  in  Table  6.  In 
these  applications,  a  match  between  observed  and  simulated  concentrations  at  source  area 
monitoring  wells  was  achieved  based  on  reasonable  estimates  of  the  mass  and  dimensions  of  the 
source  zone  along  with  composition  of  the  NAPL.  This  finding  suggests  that,  at  some  sites. 
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knowledge  of  NAPL  source  complexities  (e.g.,  interfacial  area)  is  not  required  for  a  reasonably 
accurate  estimate  of  TOR.  However,  the  dependency  of  concentration  versus  time  results  on 
source  zone  mass  and  length,  along  with  NAPL  composition,  suggests  that  source 
characterization  methods  have  the  potential  to  reduce  the  uncertainty  associated  with  source 
dissipation  calculations  using  NAS.  The  NAS  simulation  results  also  demonstrate  the  value  and 
efficiency  of  refining  TND  predictions  using  postremediation  data. 

One  finding  of  this  demonstration  was  that  NAS  proved  to  be  applicable  to  all  eight  sites, 
independent  of  hydrogeology,  contaminants,  characteristics  of  the  source  zone,  or  ERA. 
Therefore,  the  simplifying  assumptions  associated  with  the  analytical  solutions  and  the  numerical 
source  zone  model  do  not  appear  to  render  NAS  ineffective  but,  in  fact,  demonstrate  the 
applicability  and  utility  of  NAS  to  a  wide  range  of  contaminated  sites.  In  contrast, 
comprehensive  three-dimensional  numerical  models  that  are  constructed  to  simulate  the 
complexities  of  a  groundwater  system  and  features  of  a  plume  often  are  subject  to  limited  data 
and  may  include  unrealistic  boundary  conditions  that  do  not  honor  the  actual  field  conditions. 
For  example,  the  source  term  in  many  models  (e.g.,  RT3D)  is  not  based  on  the  concept  of  mass 
balance.  In  contrast,  the  source  model  in  NAS  implements  the  NAPL  dissolution  package  of 
SEAM3D  in  which  mass  is  conserved  for  all  components  in  the  aqueous,  solid,  and  NAPL 
phases.  The  results  of  this  demonstration  clearly  suggest  that  because  NAS  is  based  on  sound 
science,  it  can  serve  as  an  effective  tool  for  decision  making,  data  analysis,  and  cost  optimization 
at  a  wide  range  of  contaminated  sites  and  is  not  limited  to  a  small  subset  of  “simple  sites” 
because  of  its  simplicity.  However,  there  are  many  sites  where  complex  hydrogeology,  highly 
non-uniform  groundwater  flow,  and  the  desire  to  simulate  complicated  remediation  strategies 
will  dictate  the  use  of  a  comprehensive  numerical  model.  At  other  sites,  it  may  prove  to  be 
efficient  to  employ  NAS  as  a  precursor  to  a  comprehensive  (numerical)  solute  transport 
modeling  investigation  or  as  a  follow-up  to  a  previously  completed  modeling  study. 

6.6  END-USER  ISSUES 

Stakeholder  buy-in  is  evident  by  the  support  provided  by  Navy  Engineering  Field  Divisions.  This 
is  exemplified  by  the  recently  funded  NAS  upgrade  to  allow  for  incorporating  of  a  source 
removal  term.  Earlier  versions  of  NAS  have  been  the  subject  of  Navy-supported  conference 
presentations  (e.g.,  the  annual  RPM  Conference  in  Port  Hueneme)  and  a  Remediation  Innovative 
Technology  Seminar  (RITS)  course  module.  End-user  concerns,  reservations,  and  buy-in  factors 
all  point  to  the  RPM  willingness  to  utilize  the  NAS  approach  versus  a  more  comprehensive 
modeling  effort.  Concerns  with  implementation  have  been  addressed  through  the  NAS  site 
applications  using  feedback  from  users  during  this  demonstration. 

NAS  is  designed  for  application  to  sites  where  the  plume  is  reasonably  well-characterized, 
stable,  and  not  expanding  with  time.  NAS  is  ideal  for  application  to  sites  where  data  collection 
has  followed  DoD/EPA  protocols  for  MNA.  NAS  provides  value-added  to  sites  where  long-term 
performance  monitoring  and  redox  indicator  data  has  been  collected.  NAS  has  helped  to  address 
at  least  one  major  criticism  of  MNA  protocols.  Specifically,  technical  guidance  focused  on  what 
data  to  collect  but  did  not  adequately  explain  how  to  quantitatively  utilize  the  data  to  address 
TOR  and  other  MNA-related  remediation  questions.  NAS  provided  the  framework,  technical 
guidance,  and  computational  tool  to  achieve  this  end. 
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6.7  APPROACH  TO  REGULATORY  COMPLIANCE  AND  ACCEPTANCE 


NAS  already  is  widely  visible  to  regulatory  ageneies.  NAS  training  has  been  delivered  at  many 
venues  throughout  the  United  States  over  the  last  4  years.  This  includes  21  short  courses  to 
various  professional  organizations.  Many  have  involved  regulators,  and  in  some  cases,  the 
project  principal  investigators  have  delivered  courses  and  training  directly  to  state  and  federal 
agencies.  The  principal  investigators  participated  in  seven  invited  presentations  through  the  U.S. 
Navy  RITS  in  which  EPA  and  state  regulators  were  present.  In  general,  NAS  appears  to  be 
widely  accepted  by  regulatory  agencies  and  in  compliance  with  state  and  federal  guidelines  for 
site  remediation  investigations. 

By  documenting  the  effectiveness  of  NAS  as  a  tool  to  predict  time  of  remediation,  regulatory 
acceptance  will  be  facilitated.  Furthermore,  management  decisions  regarding  source  zone 
treatment  can  be  enhanced  through  appropriate  data  collection  activities  driven  by  model  input 
requirements  followed  by  NAS  simulations  and  cost  benefit  analyses. 
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